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Abstract 

A novel kind of irreversible phase transitions (IPT's) driven by an oscillatory 

input parameter is studied by means of computer simulations. Second order 

IPT's showing scale invariance in relevant dynamic critical properties are 

found to belong to the universality class of directed percolation. In contrast, 

the absence of universality is observed for first order IPT's. 
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Far from equilibrium system often exhibits irreversible phase transitions (IPT's) between 
an active (or reactive) regime and an inactive (or absorbing) state. Such transitions are 
irreversible because a system trapped in an absorbing state can never escape from it. Among 
others, models exhibiting IPT's are directed percolation |1J, contact processes 0, branching 
annihilating walkers ||, forest-fire models @-|6[], models for dynamic evolution of living 
individuals |§ and several models of catalyzed reactions such as the ZGB model || and 



variations 10,11 



In spite of the considerable progress achieved in the understanding of irreversible critical 



behavior (for reviews see e. g. []T0| |12j| ) , the study of the dynamic response of irreversible 
systems to external perturbations is still in its infancy. Very recently, we have studied the 
dynamic response of the ZGB model close to a second order IPT, showing that, after driving 
the system within the absorbing state, the subsequent relaxation can be well described by 
an stretched exponential decay |13| . The lack of additional studies in this field is in contrast 



to their equilibrium counterpart. In fact, the study of dynamic response systems in thermo- 
dynamic equilibrium, close to reversible phase transitions, to an external perturbation, is a 
subject of current interest [[H] . 



The aim of this work is to study the dynamic response of irreversible systems to a periodic 
external perturbation. The study has been performed close to both first and second order 
IPT's. For this purpose we have selected a forest fire model with immune trees (FFMIT) 
and the Ziff-Gulari-Barshad model (ZGB model ) ||. 

The ZGB model is a lattice gas reaction system aimed to mimic the catalytic oxidation 
of carbon monoxide, CO + (1/2)02 — > CO2, according to the Langmuir-Hinshelwood mecha- 
nism. So, reactants CO and O2 are adsorbed on the surface of the catalyst with probabilities 
Pco an d Po 2 - Since these probabilities are normalized ( Pco + Po 2 = 1)5 the ZGB model 
has a single parameter, i.e. Pco- For P C o — * 1 ( Pco - ► ), the surface of the catalyst 
becomes inactive due to complete saturation with CO ( O ) species, respectively. However, 
between this two inactive states there is a reactive regime as it is shown in Figure [I]. Close 
to P2C0 — 0.5256 the ZGB exhibits an abrupt change in the rate of reaction and reactant's 



coverages, as it is shown in Figure p], indicating a first order IPT. Further details on the 
ZGB model can be found elsewhere fll~0| . 

In order to study the dynamic response of the system, first a stationary configuration of 
the ZGB model well inside the reactive regime (actually in the center of the reaction window 
Pq Q = 0.455 ) is obtained. Subsequently, an square oscillatory perturbation of the form 

Pco + A Pco , if < t < t/2 

Pco (t) = < P™ Q , if t/2 < t < r (1) 

Pco{t) = Pco{t + t) , V t > 
is applied to the system, where Ap co is the amplitude and r is the period of the perturbation 
( see also Figure p. 

There is a great variety of forest-fire models (FFM's), e. g. for an extensive review 

see e. g. @. FFM's are stochastic cellular automata which are defined on d-dimensional 
hypercube lattices with L d sites. Each site can be either occupied by a tree, a burning tree, 
or empty. A great number of FFM's can be defined giving a set of rules which are used, 
during each time step, to update the system in parallel. These rules are: (1) burning tree — > 
empty site; (2) tree — > burning tree with probability (1 — g) if at least one nearest neighbour 
is burning; and (3) empty site — > tree with probability p. The probability p is to be taken 
as growing probability and g is the immunity of each tree to catch fire. 

In the present work we shall focus our attention to the forest-fire model with immune 
trees (FFMIT) [§-0. Qualitatively speaking, if the growing probability is non-zero (p — > 1) 
and the immunity is low (g — > 0) one expects coexistence of fire, trees and empty sites. 
However, keeping p constant and increasing g the fire will eventually cease and the system 
will become irreversibly trapped in an absorbing state with the lattice completely filled 
by trees. So, the FFMIT exhibits second order IPT's between an active state with fire 
propagation and an absorbing state where the fire becomes irreversibly extinguished ||.|6| . 

In order to study the dynamic response of the FFMIT upon temporal variations of the 
parameters, one can vary either p, g, or both of them. However, we have worked taking 
p = constant while g is varied. For this purpose, the procedure is as follows: first a stationary 



active state of the standard FFMIT is obtained for fixed values of the parameters. In this 
work we take p = p = 0.5 and g = go = 0.46. Subsequently p is kept fixed and g is varied 
harmonically according to 



where A g and T are the amplitude and the period of the oscillation, respectively. Notice 
that the critical point is given by po = Pc = 0.5 and g c = 0.5614 ± 0.0005 J|J|. 

Due to the variation of the parameters, either g(t) or Pco(t), it is expected that for 
long periods and/or large amplitudes, the systems may eventually become trapped into 
an absorbing state. So, the oscillatory variation of the parameters may cause IPT's from 
the active (oscillatory) regimes to the absorbing states. These transitions may occur at 
critical values of the amplitude (A c g and A c PcQ ) and the period (T c and r c ). In order 
to characterize and study such IPT's, we have performed epidemic simulations (ES) fD-Q. 
The idea behind ES is to start from a configuration very close to the absorbing state and 
subsequently, to follow the temporal evolution of the system under consideration. Therefore, 
we took a sample filled with trees (or CO molecules) except for a small patch of 2 x 2 sites 
having burning trees (or empty sites) placed at the center of the lattice for the FFMIT 
and the ZGB model, respectively. Depending on the values of the parameters and due to 
the stochastic nature of the processes, such a small perturbation would either propagate, or 
eventually become extinguished. During the epidemic propagation the following quantities 
are measured: (i) the average number of burning trees (empty sites) N(t), respectively; (ii) 
the survival probability P(t), i.e. the probability that the fire is still ignited (there are empty 
sites) at time t, respectively; and (iii) the average mean-square distance R 2 (t) over which 
the fire (empty sites) has spread, respectively. 

The usual ansatz for ES close to second order IPT's is to assume that N(t), P(t) and R 2 (t) 
would obey power-law dependencies with exponents r),—5 and z respectively. However, 
in the present ES of the FFMIT the input parameter varies harmonically, so we expect to 
obtain an oscillatory output modulated by a power-law, that is 

4 




(2) 



N{t)~V.co8(jrt + B), (3) 

and similarly for P(t) and R 2 {t), where B is a constant phase- shift. 

Figure 0(a) shows a log — log plot of N(t) vs. t obtained performing ES, where the 
influence of oscillatory input can clearly be observed. In order to perform a fit, we have 
first determined the values of N(t) on peaks, valleys and centers, given by N + , iV~and iV°, 
respectively. The insert of figure 0(a) shows that log — log plots of iV + , iV~and N° versus 
t can be very well fitted by straight lines with slopes r] + = 0.218 , r]~ = 0.223 , 77° = 0.224, 
respectively. As it is usual, plots drawn taking smaller (larger) amplitudes show upward 
(downward) curvature suggesting that they are off-criticality, respectively. Precisely, the 
straight line observed is the signature of a power-law behavior which characterizes a second- 
order phase transition exhibiting scale invariance. The insert in Figure Q(b) shows that it is 
also possible to determine the assumed constant phase shift in equation (|3]), comparing plots 
of N(t) and g(t) versus t. For example, in Figure 0(b), we have obtained B = 1.37 ± 0.34 . 
Therefore we can fit the whole curve of N(t) vs. t using the already determined values of rj 
and B, as it is shown in Figure 0(b). 

Figure shows a log — log plot of P(t) versus t which also exhibits oscillatory behavior. 
In this case, the drops of the survival probability are due to the increment in the immunity 
which drives the system into the absorbing state making the fire propagation harder, an 
effect which may cause the eventual extinction of some epidemics. Defining the maximum, 
medium and minimum values of P(t) in each cycle as P + , P~and P°, respectively, we can 
fit the exponent 5, as it is shown in the insert of figure 0. Our results, at criticality, are 
5 + S 0.463 , 5- S 0.400 , 5° £ 0.433 respectively. 

It is found that R 2 is less sensitive to the oscillatory input than N(t) and P(t). The plot 
of R 2 vs. t can roughly be fitted by a straight line which yields a slope z = 1.18 (not shown 
here for the sake of space). 

It should be noticed that the dynamic exponents which characterize the IPT's driven 
by the oscillatory parameter are in good agreement with those of the universality class of 
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directed percolation (DP, in 2+1 dimensions), namely 77 = 0.214, 5 = 0.460, and z = 1.134 
|Ij]. Also, the hyperscaling relation dz = 45 + 2i] [[!]], is rather well satisfied by these 
exponents. So, we conclude that the novel type of transition discussed so far, can be placed 
in the universality class of DP. This result extends the validity of Janssen's conjecture ||16|| , 



which states that a continuous transition into an absorbing state characterized by an scalar 
order parameter may belong to the universality class of DP, to second order irreversible 
transitions driven by oscillatory parameters. 

Performing ES with different values of A g and T we have evaluated the phase diagram 
of the FFMIT under oscillatory driving, as it is shown in figure f|. The critical curve T c 
vs. A 9c shows the location of second order IPT's between the active regime (trees+burning 
trees+empty sites) and the absorbing state (only trees). All of these transitions belong to the 
universality class of DP. The insert of Figure f| shows a log-log plot of T c vs. AA = A gc — (g c — 
go). Notice that AA is the "excess critical amplitude", namely a renormalized amplitude 
which accounts for the value of the oscillatory parameter which exceeds the stationary critical 
threshold allowing the system to make an excursion to the absorbing state. The data is then 
consistent with an hyperbolic-like behavior of the form T c oc AA a with exponent a = 4.55. 
The deviation from this behaviour, observed for short periods, is likely due to the fact that 
the input signal can not longer be considered as harmonic. 

Figure [5] shows log — log plots of N(t) versus t obtained performing ES of the ZGB 
model for a fixed period (r = 20 MCS ) and different values of the amplitude Ap co . In 
contrast to the FFMIT, here N(t) decreases with time as it has been observed for ES at first 
order IPT's |i~7fl . Figure [| also allows to identify subcritical (Ap co = 0.123), supercritical 



(Ap co = 0.125) and critical (Ap = 0.12345) behavior, respectively. In contrast to the ES 
of the second order IPT's of the FFMIT which exhibit power-law behavior (Figure |2[), the 
plot of figure |, at criticality, exhibit marked curvature suggesting the cross over to a cut-off 
at certain long time. Based on the observation that the short time behavior of N(t) vs. t 
seems to obey a pseudo power-law with effective exponent rfff £ 2, we have proposed |Tj| 
the following ansatz as a generalization of equation (||) 



N(t) oc (t/r*)-" e/ exp{-t/r*) , (4) 

where r* set a characteristic time scale and the oscillatory dependence of N(t) is not con- 
sidered for the sake of simplicity. Measuring the value of the average number of holes at 
half period N*(t) it is possible to check the validity of the ansatz given by equation (|]). In 
fact, the insert of Figure ^ shows a semilog arithmic plot of N(t) (t/T*) J?e vs. (t/r*) can 
be well fitted with a single parameter, namely the characteristic time r* ~ 471 ± 5. This 
results indicates that the power-law behavior for short times (t < r*) crosses over to an 
exponential decay (t > r*). Our results not only point out that the IPT's driven by the 
oscillatory parameter in the ZGB model are of first order, but also allow us to rule out the 
occurrence of power law (scale invariance) behavior. This finding conciliates the behavior of 
first order IPT's with their equilibrium counterpart where it is well know that the existence 
of short range correlations prevents the occurrence of scale invariance. 

Performing ES of the ZGB model under oscillatory driving of the input parameter for 
different values of the period we have evaluated the corresponding phase diagram, namely 
a plot of r c vs. Ap CQ , as it is shown in Figure []. The critical curve shows the location 
of first order IPT's between a reactive state with CO2 production and a poisoned inactive 
state where the catalyst's surface is fully covered by CO. The insert of Figure ^ shows a 
log-log plot of r c vs. A A = A c Pco — (P2C0 — Poo)- As i n the previous case of the FFMIT 
(Figure |j) the data is consistent with an hyperbolic-like behavior r c oc AA a with exponent 
a = 1.53. It is also interesting to know that for short periods (r c < 1 ), the value of the 
excess amplitude saturates just at AA S = Ap CQ /2 = P2C0 — Pco = 0.071. This result can 
straightforwardly be interpreted assuming that, such short periods are indeed much shorter 
than the relaxation time and consequently the system only feels the average value of the 
applied oscillatory pressure. Under these circumstances, the critical amplitude is expected 
to be Ap co = Pq + 2. AA S , as it has already been observed in Figure ||. 

In summary, IPT's induced by an oscillatory external parameter are studied in two dif- 
ferent systems, namely a forest fire model with immunity and a model for the catalytic 
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oxidation of carbon monoxide. Second order IPT's are placed in the universality class of di- 
rected percolation. However, first order IPT's lacks of universal behavior. These findings are 
in qualitative agreement with well established concepts developed in the study of reversible 
phase transitions. A phase diagram for this new type of IPT's is a critical curve, period 
versus amplitude of the oscillations which sets the boundary, between an active regime and 
an absorbing state. These kinds of phase diagrams are evaluated for both studied models 
and, in both cases, hyperbolic-like critical curves are found but the corresponding exponents 
are different. 

We expect that these findings will stimulate theoretical and experimental work on the 
field of far from equilibrium systems in general and in the study of irreversible critical 
behaviour in particular. 
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FIGURES 

FIG. 1. Plots of Rco 2 ? ®co an d @o for the ZGB model on the square lattice of side L = 512. 
In the upper part the double arrow shows the range where the CO pressure is oscillatory varied. 
More details in the text. 

FIG. 2. Results for the FFMIT obtained performing ES. (a) Log — log plot of N(t) vs. t with 
A g = 0.1695 and T = 55. The insert shows Log — log plots of N + , iV~and N° versus t (see text), 
(b) Fit of N(t) vs. t (see text). The insert shows the constant phase shift B between input and the 
output signal. 

FIG. 3. Log — log plot of P(t) versus t obtained performing ES for the FFMIT with A g = 0.1695 
and T = 55. The insert shows the linear fits giving 5 + , 5~ and 5° (for more details see the text). 

FIG. 4. Critical curve T c vs. A gc showing the location of second order IPT's driven by the oscil- 
latory input between the active regime (AR, trees+burning trees+empty sites) and the absorbing 
state (AS, only trees) in the FFMIT. The insert shows a log-log of T c versus AA = A 9c — (g c — go). 
The solid straight line has slope a = 4.55. More details in the text. 

FIG. 5. Log — log plots of N(t) versus t obtained performing ES of the ZGB model for a fixed 
period (r = 20 MCS ) and different values of Ap co as indicated in the figure. The insert shows a 
semi logarithmic plot of N(t) (t/r*) r,eff vs. t/r* . The straight line, which validates the ansatz of 
the equation @ has been fitted to the values of the output signal at half period N*(t) , obtaining 
t* = 471 ± 0.5 . More details in the text. 

FIG. 6. Critical curve r c vs. A c Pco , obtained for the ZGB model applying harmonic perturba- 
tions. The curve set the boundary between the active regime (AR) and the absorbing state (AS). 
The insert shows a log-log of r c versus AA = Ap CQ — (Pico ~ Pco)- ^ ne s °hd straight line has 
slope a = 1.527. The vertical dashed line shows the saturation of AA for r c < 1. More details in 
the text. 
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Figure 1 
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Figure 2 (a) 
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Figure 3 
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Figure 5 
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